REPORT  C 


DEFLECTIONS  OF  A  SUPERSONIC 
WING  DUE  TO  AERODYNAMIC 
HEATING 

by  Richard  W.  Luce  Jr. 

D.I.C.  Project  Number  6553 

U.S.Air  Force  Contract  Number  W33-038  oc- 17239 

JUNE  I,  1949 


AD-A278  113 


...  \ 


jJ^-ELECTEl 

P^APRaf 19941 


TW.'-  / 


94-09928 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

DIVISION  OF  INDUSTRIAL  COOPERATION 
CAMBRIDGE  •  MASSACHUSETTS 


DISCUIEI  NOTICE 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DUC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


REPORT  C 


DEFLECTIONS  OF  ▲  SDPERSONIC  WING 
DUE  TO  AERGDINAMIC  HEATING 


RICHARD  V*  LUCEy  JR* 


For  th«  U«  S«  Air  Foroo  nnder 
Contract  W33-038  aool7239 
D.I.C.  Project  No* 


Jane  1,  1949 


DinSKNI  OF  INDDSTRIAL  COOPERATION 
MASSACHUSETTS  INSTITUTE  OP  TECHNOLOGI 
CAMBRIDGE,  MASSACHUSETTS 


zmm 


This  report  is  one  of  a  series  irtiieh  is  the  result  of  a  joint 
undertaking  of  the  Depcurtaents  of  Aeronautical  Engineering,  under 
the  8upenri8i<»i  of  Professors  J.  S.  Newell  and  H.  S.  Tsien,  and  of 
Mechanical  Engineering,  under  the  supervision  of  Professors  J.  Kaye 
and  W.  M.  Rohsenow,  at  the  Massat^usetts  Institute  of  Technology  for 
the  U,  S.  Air  Forces  under  Contract  No.  1733*038  ac*17239*  The  work 
was  perfozmed  in  the  Division  of  Industrial  Cooperation  as  projects 
Die  6553  and  DIC  6580.  The  work  began  in  October,  1947  and  extended 
to  June,  1949* 

Since  the  basic  problems  relating  to  the  thermal  effects  in  a 
structure  moving  at  supersonic  speeds  are  relatively  new  and  unsolved, 
the  major  efforts  in  this  study  were  of  an  exploratory  nature.  Rather 
than  attempting  to  cover  a  large  number  of  different  structures,  it 
was  deemed  advisable  to  select  a  particular  case,  namely  that  of  a 
thin  wedge*shape  wing  travelling  at  supersonic  speeds,  and  to  dis¬ 
cover  the  basic  problems  which  must  be  solved  in  order  to  design  such 
a  structure  safely.  Even  for  this  particular  case  only  the  more  im¬ 
portant  variables  could  be  examined  in  the  time  available. 

The  problem  was  divided  into  the  following  four  parts: 

a.  The  determination  of  the  temperature  distribution  in  the 
wedge-shape  wing.  The  results  are  given  in  Report  A. 

b.  The  thermal  stresses  in  the  wing  caused  by  this  non- 
tiniform  distribution.  The  results  are  given  in  Report  B. 

c.  The  effects  of  this  temperature  distribution  on  the  aero¬ 
dynamic  behavior  of  the  wing.  The  results  are  given  in 
Report  C. 

d.  A  search  of  the  literature  in  order  to  discover  any  previous 
.work  in  this  field.  The  results  are  given  in  Report  D. 


In  thla  report,  a  stody  la  nada  of  the  defleetlona  of  a 


wing  due  to  aerodynanlo  heating,  and  the  effeeta  of  theae  defleetlona 
on  the  aerodynanlo  perfomanee  of  the  wing.  The  te^>erature  diatribu- 
ti(»  used  la  that  calculated  In  Report  A  of  thla  seriea  hy  a  Method 
of  finite  dlfferencea.  A  different  approach  to  the  heat  tranafer 
problea  la  auggeated  in  an  Appendix  to  this  Report.  For  the  sake  of 
aiaplleity,  the  wing  whldi  is  considered  is  a  flat  rectangular  plate 
of  constant  thickness.  Rowewer,  methods  of  generalising  the  procedure 
to  solid  wings  of  Tariable  thickness  are  indicated. 

The  analysis  is  baaed  upon  a  strain  energy  procedure  rather 
than  a  solution  of  the  differential  equations  of  themo-elasticlty. 

In  this  procedure,  it  was  found  necessary  to  take  into  account  the 
coupling  of  stresses  in  the  plane  of  the  plate  with  transrerse  deflec¬ 
tions  of  ^e  plate. 

Results  are  given  in  tabular  and  graphical  fom,  and  show 
that  in  general  the  effect  of  these  defleotims  is  not  large.  The 
principal  aerodynamic  effect  is  an  increase  in  drag  coefficient.  How¬ 
ever  all  of  these  calculations  are  based  upon  a  llach  number  of  six. 
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As  shown  by  ths  results  of  calculations  in  the  previous  reports 
of  this  seriesf  non-uniform  heating  by  rapid  acceleration  to  supersonic 
speeds  in  dense  atmosphere  leads  to  very  large  thermal  stresses  in  the 
material.  These  thermal  stresses  are  calculated  on  the  assumption  that 
the  dimensions  of  the  wing  are  infinitely  large.  If  the  dimensions  of 
the  wing  are  finite,  there  will  be  some  release  of  stresses  due  to  the 
effects  of  free  edges.  To  determine  this  effect  is  part  of  the  aim  of 
this  report.  However,  since  the  material  will  be  assumed  to  be  con^ltte^y 
elastic,  idiile  plasticity  of  the  material  must  necessarily  appear  at  the 
high  prevailing  temperature,  the  magnitude  of  the  stresses  so  calculated 
cannot  be  considered  as  quantitative.  The  main  purpose  of  the  present 
study  is  rather  the  determination  of  the  deformation  of  the  wing  due  to 
aerodynamic  heating.  As  will  be  seen  later,  the  deformation  is  indepen'- 
dent  of  the  value  of  Young's  modulus  and  is  only  a  function  of  the  tem¬ 
perature  distribution  and  the  thermal  coefficient  of  expansion  of  the 
material.  Therefore,  the  appearance  of  plasticity  at  high  temperature 
and  the  resultant  changes  in  the  effective  elastic  modulus,  will  not 
essentially  modify  the  calculated  results  on  deformation.  Hence,  the 
assumption  of  complete  elasticity  is  justified  for  the  calculation  of 
the  deformation  of  the  wing. 

Specifically  then,  this  report  attempts  to  obtain  a  first 
order  estimate  of  the  transverse  deflection  of  a  cantilever  wing  caused 
by  the  thermal  stresses  set  up  by  a  transient  teB;>erature  distribution. 

The  wing  is  assumed  to  be  accelerated  at  a  constant  rate  starting  at  a 
llach  number  of  1.4.  The  acceleration  is  assumed  to  occur  at  constant 
altitude  and  constant  wing  loading.  The  angle  pf  attack  is  thus  deter¬ 
mined  to  maintain  this  constant  wing  loading  during  the  period  of  accel¬ 
eration.  During  this  period  of  acceleration  there  is  heat  transfer  froa 
the  air  into  the  wing,  producing  a  tenqperature  distribution  as  a  function 
of  time.  This  temperat\ire  distribution  has  been  calculated  for  a  partic¬ 
ular  set  of  conditions  by  a  method  of  finite  differences  as  given  in  Report 
A  of  this  series  and  it  is  this  temperature  distribution  which  will  be  used 
in  the  following  work.  ♦ 

To  calculate  the  deformation  and  the  stresses  with  a  given  tem¬ 
perature  distribution,  one  can  proceed  by  either  of  two  ways:  1)  To 
solve  the  exact  equilibrium  equations  and  the  equation  of  canpatiblllty; 

2)  To  use  the  energy  method.  This  first  method  is  more  exact  but  the 
labor  Involved  in  actual  computation  is  prohibitive  for  the  case  of  a 
cantilever  wing.  The  second  method  will  thus  be  used.  This  method  has 
the  advantage  of  being  flexible  and  can  be  adapted  without  any  difficulty 
to  wings  with  arbltrau^  thickness  variation.  For  the  detailed  analysis 
given  below,  the  thickness  of  the  wing  is  assumed  to  be  constant  to  sim¬ 
plify  the  calculation.  This  is  believed  to  be  justified  for  t*e  present 
study,  as  the  aim  is  to  obtain  a  first  order  estimate  of  the  deformation 
rather  than  the  exact  numerical  value. 


*In  an  Appendix  to  the  present  Report,  an  alternate  method  based  upon  the 
Laplace  Transform  will  be  presented. 
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NOTATIORS 

1)  ■  local  drag  coefficient 

2)  Cj[  ■  local  lift  coefficient 

3)  ■  local  moment  coefficient 

4)  s  8  pressure  coefficient 

5)  "D  :  *  plate  bending  stiflbess 

6)  £  8  Young's  Modulus  of  Elasticity 

S  =  -zlxv) 

d)  K  s  plate  thickness 

9)  j2.  s  semichord 

10)  L.  s  semispan 

11)  ti  )  z  Laplace  transform  of  (  } 

12) N,»v»oi  z  force  in  ac-direction  per  foot 

13)  force  in  y-direction  per  foot 

14)  Na^  *  H  *  shear  load  per  foot 

15)  R  »  Vjt 

16)  S  •  Vl. 

17)  T  m  temperature 

IS)  *T^  ■  temperature  averaged  with  respect  of  Z 

19)  H.  s  displacement  in  the  x~direotion 

20)  V*  s  di^lacement  in  the  y-directioc 

21)  V  ■  potential  energy 

22)  ur  s  displacement  in  the  Z  directloo 

23)  W  a  strain  energy  function 

24)  %  ■  spanwlse  coordinate 

25)  8  chordwise  coordinate 


26}  Z  s  transverse  coordinate 

27)  OC  r  linear  thermal  coefficient  of  expansion 

28)  shearing  strain 

29)  Eji  r  tensile  strain  In  the  x-dlrectlon 

30)  ■  tensile  strain  In  the  y-direction 

31)  Til  «  Va 

32)  r  time 

33)  V  s  Poisson’s  ratio 

34)  \  I  */i, 

35)  TT  s  strain  energy  of  a  flat  plate 

36)  TTi  >  strain  energy  due  to  bending  stresses 

37)  TTx  »  strain  energy  due  to  plane  stresses 

38)  Oj^  s  tensile  stress  In  the  x-dlreetlcm 

39)  r  tensile  stress  In  the  y-direction 

40)  shearing  stress 

41)  ^  s  thickness  ratio  ■  Vm 

42)  ^  r  Airy's  stress  function 

FORMULATION  OF  THE  PROBLEM 


The  coordinate  sj^stem  and  dimensions  to  be  used  in  the  analysis  are 
as  diomn  in  Figure  1.  The  followiag  set  of  non-dimensional  parameters  will 
also  be  used  in  the  analysis. 


The  results  of  the  finite  difference  procedure  indicate*  ttiat.  the 
temperature  distribution  at  any  instant  can  be  approximated  quite  closely  by 
an  expression  of  the  following  form,  where  ^  i ’Z)  Is  the  local  teo^er- 

ature  at  the  point  (X,  ^,^)  above  the  uniform  temperature  before  heating. 


TCX.»,Z1- AW-*-  -►  U) 


*See  Report  C 
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The  eoefflolents  A,  and  C  are  only  functions  of  the  span-vise  coordinate  y. 
The  influence  of  the  wing  tips  on  the  teo^rature  distribution  is  thus  neglected. 

Frcn  this  espressionf  it  is  easily  seen  that  A,  B,  and  C  have  the  fol- 
lovisg  significance. 


Ai>«  =  t;., 

Bw- 

It  is  also  found  convenient  to  define  the  mean  temperature  as: 

=  “K  \  «  AUW 

^  ^  A 


The  deflections  of  a  wing  are  mainly  due  to  the  bending  of  the  wing 
surface.  But  bending  is  the  result  of  the  difference  in  e3q}an8ion  or  the 
difference  in  temperature  of  the  material  at  the  top  and  the  bottom  svirfaces. 
If  this  is  true,  then  only  the  linear  term  ^  temperature  dis¬ 

tribution  vlll  influence  the  deflectims.  However,  this  is  not  true:  Due  to 
the  non-uniformity  of  the  mean  temperature  defined  in  (4);  there  is  rather 
large  stretching  and  con^ression  of  the  median  surface,  lying  half-way  between 
the  top  and  the  bottom  surfaces.  The  stretdiing  and  ccxupresslcm  of  the  median 
surface  yield  an  average  tensile  or  e(m9>re8sive  stress  in  the  wing.  These  are 
the  stressed  in  the  plane  of  the  wing.  Such  stresses  will  Influence  the  bend¬ 
ing  of  the  wing  through  coupling  effects.  Therefore,  to  calculate  the  deflec¬ 
tion  of  the  wing,  the  procedure  can  be  divided  into  two  steps: 

1)  Calculate  the  average  or  plate  stresses 

2)  Calculate  the  bending  with  the  coupling  effects  of  the  plate 
stresses. 

T^ese  steps  will  be  treated  separately  in  the  following  sections. 

PLATE  STRE&J:S 


Let 


r  tensile  strain  in  the  x-direction 
£y=  tensile  strain  in  the  y-dlrection 
'lx,-  shearing  strain 
QTii  s  tensile  stress  in  the  x-directlon 
O^s  tensile  stress  in  the  y-dlrection 


5 


sheariiig  stress 
C  s  Young's  nodulus 
y  5  Poisson's  ratio 

OC  s  linear  themal  expansion  coefficient. 
Then  the  relations  between  the  strain  and  stress  are: 

€.  -  ■  '’  °»)  +  « 

h  =  e"^“»  “  ''*'»)  * 

a  _  •r- 

0«JJ - g— t 


The  equilibrium  equations  are: 


ISl 


In  addition,  the  cooq}atibility  equations  becone: 


16) 


(7) 


Combining  the  stress-strain  equations  and  the  compatibility  equation: 

^(ok-va-,)  H.  ^(r, - voi)  +  Esc  V*t;.  =  an+v)  C») 


But  differentiating  the  equilibrium  equations  mith  respect  to  x  and  y  respec¬ 
tively,  and  adding: 


— 


(9) 
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H«nce  toy  eliminatlcg 

Vft-v  ^  ^  ^  ^  ^  Etc  V‘X. 


or: 


^‘Coi  ■>-<ii)=-EolV'-TUv 

Now  introduce  the  stress  function  defined  toy: 

Ox 

1)  becomes: 

^  *  -  E  cC  v*-  TU. 


lit 

Tv^  1 


!*• 

w  ^ 


Tx^  = 


i 


(\6) 


OO 


02) 


t\3) 


Thus  the  differential  equation  (V^)  must  be  solved  subject  to  the 
boundary  conditions  of  zero  normal  and  shear  stress  on  the  free  edges.  The 
supported  end  of  the  plate  is  assumed  to  be  free  to  displace  in  the  y-diree- 
tlon.  Tha  same  boundary  conditions  are  satisfied  if  the  cantilever  plate  is 
considered  as  one  half  of  a  rectangular  plate  free  at  all  edges.  Thus  for 
convenience,  the  cantilever  plate  is  replaced  by  an  equivalent  free  plate  of 
twice  the  length,  as  shown  in  Figure  2. 

The  boundary  conditions  became: 


Ae  *  o 

A't  X  »  iU.  j  ■»  «  O 


A  particular  solution  to  equation  (l)  is: 

<t>p  = -EK  A-iAa-  UJ) 

^  S-  O 

Tv-i.. 


Then: 
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The  problem  Is  thus  reduced  to  that  of  finding  a  homogeneous  solu¬ 
tion  to  the  equation: 


With  the  boundary  conditions: 

At  X  «  *  L.  ^  O 

-ecCTU  07) 

At  ^ai  ^HXf  *  '*■  ® 

But  this  system  is  identical  with  that  of  finding  the  stress  dis¬ 
tribution  in  a  plate  loaded  as  shown  in  Figure  3  • 

SOLUTION  OF  THE  FLAME  STRF:R.q  PROBLEM 

This  problem  can  be  solved  by  a  strain  energy  procedure  developed 
for  plane  stress.  The  problem  is  equivalent  to  that  of  finding  a  function 
which  satisfies  the  boundary  conditions  and  minimizes  the  integral:  * 

^  ®  IE  (♦niMii*’  *+  AX  (I  B) 

A  solution  is  assiuned  in  the  form. 


where  the  0(  ^  are  undetermined  coefficients. 

It  can  easily  be  seen  that  this  function  automatically  satisfies  the 
boundary  conditions)  and  hence  it  is  only  necessary  to  minimize  the  strain 
energy.  It  is  assumed  that  the  mean  temperature  can  be  expressed  in  a  power 
series  as  given  below: 

-  II 

x»e  ' 


» 


♦See  for  instance:  Timoshenko,  "Theory  of  Elasticity",  McGraw-Hill,  pp.  151 


Substituting  these  expressicms  into  the  expression  for  strain  enersy. 
there  follows: 

■Jf  EV  a 

•¥-\\o.A>V^  jr.  SStl*-  lo.^:QoC,^  -V  C*-*’*^^  -►.S'eiR*' 

-4.  ■♦•  J.3C2  R'  -*•  .  896]  ocj 


where  R  is  the  "aspect  ratio" 

But  the  minimizing  conditions  are  that: 

,o 


Wll-*  »e.«H*-|.a!6  8R*-  •♦.|.f«*|  So 

[13. r  R* +  11.11  +{l.'J73R*  + 

— R*  HI '  ♦  -  ■^) 

■Hie 

3.775] ^a.  -4-][z.«43'  R*^  -4-  6.5“ 7  -4- T5.  ij 


U(a 


C2-3b) 


'  BTamDla.  the  wing  frtiose  temperature  distribution  was  calcu¬ 

lated  before  by  the  finite  difference  method,*  will  be  analyzed.  This  wing 
has  a  thickness  h  of  0.22  ft.  and  a  semi-chord  JL  of  2.75  ft.  The  thick¬ 
ness  ratio  is  thus  8%,  The  wing  is  loaded  to  100  lbs.  per  square  foot  and 
is  accelerated  at  50,000  ft.  altitude  from  Uach  number  1.4  to  liach  ntimber  6 
at  1  Mach  number  every  30  seconds.  The  calculated  temperatur-  distribution 
at  ISaoh  number  6  is  given  in  Report  A.  Eicpressing  this  temperature  distribu 
tion  as  a  polynomial  as  expressed  by  (C>o)  $  the  following  coefficients  are 
determined ; 


•t.  >  317 
•t,  s.m.2 
252 
•bj  =  765 
2135 
=  2638 
ta  s  0 
s  -2057 

Using  a  value  of  R  «  3,  the  three  equations  becc^e: 


uz) 


IBI3 


>3*4.  9 


St. 


V5-4.3 

5Z.\ 


43*,  aroo 

~  65*00 


whose  solutions  are: 


OCa 

CSCl?- 


-  3.  9.7L 

-  js.r 

-65:4- 


U4) 


*In  an  Appendix  to  the  present  Report,  an  alternate  method  based  upon  the 
Laplace  transform  will  be  presented. 
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Using  those  valuesi  the  stress  function  beeches: 

aecdJo'-O'Cn'-O't'S  ®-  e5.4's‘-  a.zaiil  us-^ 


The  stresses  then  are, 


ff  ^**4  4_  4*4. 

a  EoC R*(  SlCl -'X'»*'-  0’’  I- 3.ZZ) 

+  •K.V-iXV-'f 


TT,.  »_  3^4. 

•»X»J. 


-E«R^-3.aa(4'iX’i'-i)(_T^f-_i-^*  +  4Tnln‘-iXV’-0'<r''»o.8\') 

4  68^41*--  3.za^)l  tZ8) 


The  results  of  the  calculation  are  presented  In  Fig.ff  ,  It  is  seen 
that  only  the  direct  stress  In  the  span-wise  direction  or  Oi  is  large.  This 
is  to  be  expected  as  the  large  temperature  changes  in  the  chord-wise  direction 
would  require  large  changes  in  the  span-wise  expansion  along  a  chord.  But  such 
large  variations  are  prevented  by  the  continuity  of  the  material  in  the  span- 
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vise  direction.  This  explains  large  values  for  the  span-wise  direct  stress. 
Kear  the  leading  edge  and  the  trailing  edge,  where  the  temperature  is  high, 
the  material  tends  to  expand  more  than  allowed,  thus  Ok  is  negative  or  com¬ 
pressive.  Kear  the  mid-chord,  the  temperature  is  low;  thus  the  material  tends 
to  expand  less  than  allowed,  the  O^is  positive  or  the  stress  is  a  tension. 

Of  course,  at  the  wing  tip,  ^  *  1,  the  direct  stress  has  to  be  zero;  there¬ 
fore  Ok  must  decrease  as  Increases  from  ^  ■  0  to\  ■  1.  This  decrease  is 
however  rather  slow  for  small  values  of  5  •  Hence,  the  tip  effect  is  limited 
to  a  rather  small  region. 

It  should  be  noted  that  the  above  result  is  obtained  without  any 
assxunption  as  to  the  values  fo  £  and  OC  .  In  fact  the  stresses  are  directly 
proportional  to  the  product  E  9.  .  The  only  other  explicit  parameter  is  the 
aspect  ratio  R  of  the  wing.  However,  one  must  bear  in  mind  the  effect  of 
scale  on  the  heat  transfer  and  thus  the  calculated  result  is  not  strictly 
applicable  to  wings  of  other  dimensions.  The  influence  of  R  can  be  easily 
caiculated  by  going  back  to  (it) . 

BEKDING  DEFLECTION 

Having  calculated  t'^e  plane  stress  distribution,  it  is  now  desired 
to  introduce  it  into  the  calculation  of  the  transverse  deflections  of  the 
plate.  To  do  this,  the  total  strain  energy  stored  in  the  plate  as  it  deflects 
is  divided  into  two  parts;  first  that  due  to  the  components  of  thermal  stress 
producing  pure  bending  moments;  and  second  that  due  to  the  coupling  of  tran¬ 
sverse  deflections  and  plane  stress.  It  is  the  sum  of  these  two  energy  expo¬ 
nents  irtiich  is  then  to  be  minimized.  First  the  energy  due  to  the  bending 
stresses  will  be  developed. 

To  calculate  this  expression,  the  usual  plate  assumptions  on  the 
stresses  are  made.  These  assumptions  are  that: 


9 


4C  o; ,  o;  ,  T,, 


Accordingly,  the  equilibrium  equations  and  stress-strain  equations 

beexe: 
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The  first,  second  and  fourth  of  the  stress*strain  equations  can  be 
reuritten  as: 


Oi  -f  y  o-^v^o*.-f| 

-vv  £x-<.i^y')«^T3 
*  G 


The  next  assumption  is  that  the  deflection  ur  may  be  approximated 
by  ,  the  deflection  of  the  mid  plane  of  the  plate.  With  this 

assumption,  the  last  two  stress  strain  equations  can  be  integrated  directly 
to  yield: 


tC  a 


V  a  — 


)«• 


(831 


It  can  be  seen  that  these  expressions  are  equivalent  to  the  usual 
bending  assumptions  of  displacement  proportional  to  distance  from  the  neutral 
axis. 
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Thus  integrating  this  function  directly  with  respect  to  Z: 

where  B  is  the  coefficient  of  the  linear  term  in  ( 2 ) . 

The  expression  for  the  strain  energy  due  to  the  coupling  of  plane 
stress  and  transverse  deflections  is:* 

Nhere;  m  K  0^ 

Nj  ■  Kaj 

V» 

The  total  energy  of  the  plate  is  thus  the  sum  of  these  two  compo¬ 
nents  >  namely: 

Tf  -  TT,  +  Tf  X 

The  eq''.lllbrlum  position  of  the  plate  is  determined  by  minimizing  the  energy 

tr  . 


♦See  Timoshenko  in  his  "Theory  of  Plates  and  Shells",  pp  305 
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SOLUTIOH  FOR  THE  BOiDIKG  DEFLECTIOH 

It  is  now  necessary  to  assuse  a  solution  for  the  diqplaeeneat 
shieh  satisfies  the  necessary  boundary  conditions  <m  the  edges  idiere 
dlsplacenents  are  prescribed,  and  vhidi  ccmtains  \uideterained  coeffi¬ 
cients  to  be  determined  by  the  rariational  procedure.  In  the  ease  of 
the  cantilever  plate,  the  only  edge  on  «hl(^  dlsplacenents  are  pre¬ 
scribed  is  the  clasqped  root,  and  hence  the  free  edge  boundary  conditions 
do  not  need  to  be  eaqpllcitly  satisfied  in  the  fom  of  the  assnned  deflec¬ 
tion.  This  faot,  idiloh  can  be  substantiated  in  Treffts  for  example*,  does 
not  seen  to  be  generally  mentioned  in  well-knonn  texts  on  the  theory  of 
elasticity.  Physically  this  is  easy  to  understand,  as  a  free  edge  would 
certainly  seek  the  position  of  ninlmun  strain  energy  among  all  klnena- 
tlcally  possible  forms  of  deflection.  This  situation  is  of  fundamental 
importance,  because  it  is  the  satisfying  of  the  free  edge  boundary  condi¬ 
tions  whloh  makes  the  solutlcm  of  the  cantilever  plate  problem  difficult 
to  obtain  by  other  methods.  A  form  of  the  deflection  which  meets  these 
requirements,  and  which  will  be  used  in  the  variational  procedure,  is  as 
follows: 


Then  the  second  group  of  terms  represents  the  diange  in  angle  of  attack  of 
the  section  and  the  third  group  of  te'ms  represents  the  change  in  the  curva¬ 
ture  of  the  section.  Introducing  this  expression,  and  the  previously 
calculated  Ott  ^  ,  Anp  into  the  strain  energy  expression,  and 

assuming  that: 

%  =  £  «•*■»(■*' 

4|aO 


where  the  have  to  be  determined  from  the  given  temperature  distribution: 

+  0,f  b|  b|  ^  ^  (L|  OLj  8 


•Treffts:  llathematische  Elastiaitatstheorie .  Handbuch  der  Physik, 
Tol.  6,  Kap.  2,  pp.  66,  (1926) 
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4.  4.Sa»W,  *'V^ 

■•■^  —  H4>i4,(  —  ff  Wi*-*  —  a«.i  b» 

-^b|bi3  ■*•  ■*■9  b.) 

A»*  ^ 

*  (SA  +  3  ^  3  ^  R»A)  W,]"^ 

■’■•■***A«.')  •►SbitlA,  -4.  O.IA*,*) 

-V  (  \  .©Jtl  A|  —  .  11  Ak,') 

4-a.,b,t-\.o4*?  A|-.i4-ia  A*.')  ^  a-jWiC**  l  A, 

-V  b,b,(,'7-4sA.-.ia5S’A»)  ■►«.^(l.3oA,  -*• .  1 1  Cfi- A  a.) 
+  *'*b^^\.o4’i A|+ .♦ilAa)  •vb^^\.»*^Ai  •v.ini  A*,) 


Ai-.oai  Aa)  4.9  bs^-  l.iK  A,-  .UliAa^ 

4-  bK'  Aflt  A|  H.  .  168  AO  ba.  t-.o8i  Aa) 


+  bn. bj  .  o 3 8  A3 
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froB  the  plane  streee  calculations  . 

As  a  concrete  exanole.  the  wing  considered  In  the  Section  on 
plane  stress  will  be  used  again.  Froa  previous  tenperature  calculationsy 
the  coefficients  in  the  expansion  for  ere  as  follows: 


OC,  -a  —  X.4.0 

0C4.  -a. 

OA,  »  —  3*1.6 

0C3.  «  35‘3".  6 

«i  a  —  3‘54- 

OCg  m  l3-®3 

Cig  a  —  a04- 

Assuning  the  wing  is  Bade  of  steel.  Introducing  a  value  of  oCa  fi. 
y«.S^and  a  value  of  R  s  3»  the  equation  becones. 


S*-  -v  »a.l3a,a.,  S* 

4-40.-V  4  It.a  +  4Z  w,  s» 

•»>34.acL\s^  +  aa-fi  4.a.bs,  S**  4-  3-3.  J  bj  s* 

4- 34.-4  «».S  S*-  4*l3L3a,b,S‘ 

—  .04a.5'3-#i  (xc,  4-3b|)  —  .  ••e/S4A  ♦  JA*) 

—  .  0/7/Z  6fX97  b 

+  ^-5’0.3aj*-  31.3X0.,  b,  — -^.r.oX  b*-  -  5*.  73  a.,*. 


•baalA.,  bg  4*  ^8.70  0.3  b,  4* 


.  8o  b I  bg 
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—  33.^^  —  -^S.70  o-x  b*  —  “i*  0.7  8  b;^ 

«-Z«.78ai  7-  4A.9Z  aj  b»  —  3-^.  17 

4-  1.775*  a,  4^4  -►  «..  1ST  bx 


(46) 


Coiibining  terms j  this  beeomeSi 

5^  =(ia,3«0  3’'- J0.3)a.t  +  (38,300  S*- 31. 32)*,  b, 
+  (38,900  S’-—  +5-.02)l^-  +-(31,  9ffO  S*-  5.73)  a, 
+(65,500  S’- -4-  +«.7o)b,a,^  +(30j+3a  S’  3e.i)  b,A, 
+.(S8.000S*  +  +.eo)  b,ba  +(56,750  S’-- 39. 5;^o.^ 
+(i3r,*-.oa’--  +».7.)a..  b,.  +(86,soo  8*-4«.7s)  b'^ 

+  (l53,800S’-2.8.78)a‘  +(310,r..9’-  +4X.8z)a,b, 
+  (3*3,0003’-  3-+.17)  Wi  +(t.77j-)  a.,b, 

+  (3.l55)b»b»  — +.lt7m2.q,,  +  3  b,) 

-  ,3175  h,  (3«,»  +  3  b*)  —  Z7. 9  A  a,  _  Zj.o  k  6, 


(A-n) 
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Dlfferentlatixig  this  e39)re88ian  with  re8peet  to  the  eoefflelent8  gives  the  8iz 
equations  to  detemine  the  Amo  Ci  : 

5**-  SL3Z)  b»  -v  C»l,85-oS*- J.7s)«i 
4.(ao,A>80d*‘ H-ao.Obj  «  %.ZS‘AW 

^aljSooS*-  3l.3Z^a,  •»-(77,800S*‘-.4.3-.oa)b,  +t^5.5-ooS*4.4-g.7)tt.3 

**■ ‘^•8o)b3  «  IZ.SSK. 

(in,3-«oS*-73.O-0s.  -^-OBr^roo  S*-  »  .  «3S  k 

<iir,jroo  s*-- 4.1.7)«.4.  -►^l7Z,CooS‘-«l.5X)b*  7750-3 

*4-  Z.I55  b3  *.  35Z5  k 

(8l,350S‘-5'.73)a,  •►CC5,5g«S* -►-♦t.7)b,  7-  1.775  b,. 
*^(lo€,oooS*-57.  5C)q*3  •vt3IO,5ooS''^4'Z.3x)b3  *5t  Z7.9  k 

(30,430 3*’ -v  SO. a,  -*-(63,000  5*  +  ♦.So)  b,  3.156  b*. 

4- (31 0,600  5*  4  4 3. 9 3)0*3  4- 6, 000  5*-  5  6.  3  4-)  «■  Z1 .0  k 


As  jet;  no  asstuqptions  on  the  geometry  of  the  plate  have  been  made 
other  than  R  :  3.  The  e3q>re88ion8  for  the  ten^erature  distribution,  however, 
were  derived  for  a  particular  chord  and  thiobiess  ratio,  and  presunably  would 
be  different  for  other  sections*  The  values  used  in  the  temperature  calcula¬ 
tions  were: 
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1  »  Z.TS  F-n 


Since  R  has  been  q[}eeifled  as  3, then  L  s  S.25  ft*  and  hence  S  a  .0267.  For 
these  values^  the  coefficients  satisfying  tba  six  equations  are: 


a,  «  o^o^3 

b,  s  .  0]£« 

aa,m  .oolim  (4.9) 

b^a  .  00X738 
a^m  .  oottf 
b,«  .ooBsra 

Hence,  the  deflection  ur  in  feet  beocaies: 


ur 


[.OOU8\*-*-.OOXTSt’l*]‘>| 


4-  00  s  1  r  4- .  oo  8rs  >1*' 
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A  plot  of  this  calculated  deflection  is  given  in  Figure  6.  It  is 
seen  that  the  deflection  is  greatest  at  the  wing  tip.  However,  even  at  the 
wing  tip  the  deflection  is  quite  snail.  The  significance  of  these  results 
will  be  discussed  in  a  later  Section* 

METHOD  OF  CALCDLATION  FOR  SOLID  WINGS  OF  VARIABLE  THTCKNTgSS 

The  procedure  developed  in  this  work  is  applicable  to  wings  of 
variable  thickness  provided  that  several  aodlfieations  are  nade.  The  only 
(dwnges  involve  different  limits  of  integration  in  the  expressions  for  the 
various  strain  energies.  In  all  eases,  the  thickness  oust  be  left 

inside  the  integral  when  integrating  with  reqiect  to  x,  and  7.  The  aodiflea* 
tions  are  as  follows: 


Equation  (18)  beeonesi 


Equation  (38)  boeoneat 

-eaCi-W)(urK*a  - 

equation  (39)  beeoneet 

KtNtt  {o;  ( W  +  <>»iV?r 

+*-'''»»tV5H^)}Ai'«»-»  C39*.) 

Using  these  nodified  expressians*  the  rest  of  the  procedure  aay  be 
carried  cut  exactly  as  in  the  case  of  the  constant  thickness  wing. 

GEHERAL  CONCLUSIOKS  ON  THE  DEFLECTION  BY  AEKODYBAMIC  HEiTING 

For  a  ulng  of  thickness  equal  to  0*22  ft.,  chord  2.75  ft.,  accele¬ 
rated  at  an  altitude  of  50,000  ft.,  100  lbs  per  ft.^  loading  ahd  1  Mach  nuaber 
erery  30  seconds,  the  deflections  due  to  aerodynamic  heating  at  Mach  nuaber  6 
are  given  in  Table  1  for  various  values  of  the  coefficient  of  linear  theraal 
expansions  and  aspect  ratios.  The  aethod  used  in  making  these  calcu^tions  is 
explained  in  the  previous  Sections.  This  Table  also  gives  the  changes  in  the 
aerodynaaic  coefficients  of  the  wing  due  to  such  deflections.  The  aerodynaaic 
coefficients  are  calculated  by  using  the  simple  Ackeret  formulae*  for  thin  wings 
It  is  seen  that  in  all  eases  the  change  in  the  aerodynamic  coefficient  are  saall 

Therefore,  the  results  of  the  present  analysis'  indicate  that  the 
deflection  and  the  effects  of  deflections  due  to  aerodynamic  heatixig  on  the 
aerodynamic  performance  will,  in  general,  be  small. 

This  conclusion  is  further  strengthened  by  the  following  observations: 

1)  The  calculated  deflections  and  the  effects  of  deflections  are  Inde¬ 
pendent  of  the  Young's  modulus.  Therefore,  the  decrease  of  Young's 
modulus  at  high  temperature  of  the  material  will  not  influence  the 
results  of  ealeulati<». 


*  See  Appendix  II 
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2)  At  hinder  tao^eratures,  the  Young's  nodulus  will  bs  deoraased. 
Than  aith  a  givan  aarodynande  loading,  the  daflaetiona  due  to 
this  aarodynanlo  load  will  ba  increased  at  higher  taaperaturas. 
Therefore,  at  the  hl{^  tenparature  caused  by  aerodynamic  heating, 
the  deflections  due  to  the  aarodynanlc  load  will  be  much  larger 
than  the  deflections  due  to  non’^ifozm  heating. 

3)  If  the  aircraft  is  flown  at  a  mudi  hi^^er  altitude,  the  lower 
air  densi^  will  greatly  reduce  the  rate  of  heat  input  to  the 
wing.  Then  the  non^unlformity  of  temperature  distribution  orer 
the  material  of  the  wing  will  be  greatly  reduced.  This  will  in 
turn  greatly  decrease  the  importance  of  the  deflections  due  to 
aerodynamic  heating. 
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A  fdng  is  assuMd  to  aoeslerata  at  a  constant  rata  starting  froa 
steady  flight  at  soae  super sonic  liadh  number.  The  adiabatic  mall  tempera¬ 
ture  mhieh  acts  on  the  surface  of  the  ming  is  thus  a  ftmctlon  of  time,  and 
a  non-steady  state  heat  transfer  problem  arises.  It  is  assumed  that  there 
is  no  heat  flow  in  the  spanmise  direction  of  the  ming,  and  thus  the  problem 
is  a  tmo  dlmenslcoal  one.  The  coordinate  system  used  on  the  section  is  as 
shomn  in  Figure  4* 


The  differential  equation  of  heat  transfer  for  this  system  is: 


.  ^*1“  _  I  VT 

•hT.**  “  ot 


Since  ^  ^y*and  ^''T^m'are  of  the  same  order  of  magnitude, 
for  small  ralues  of  *tr  :  ' 


4-^ 


%  VlT 


i  <r 


Hence,  assuming  that 

ential  equation  becomes:  ' 


is  negligible,  the  dlffer- 


Vr  XL 

"hz*-  >• 
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lAiioh  la  the  ona-diMnalonal  equation. 

Boundary  ecmdltlona:  The  plate  la  aaaumed  to  be  Inaulated  at  the  aorfaee 
Z  a  0,  and  heated  according  to  Newton 'a  law  at  the  aorfaee  tz.  a  k  •  At 
tlM^s  0,  the  teaperature  la  aaauaed  to  be  Identically  aero.  Ttaa  the 
boundary  condltlona  are: 

-e-  *  o  *  T  «  o 
z«o  t 

^  ^  (TV-r  —  T) 

It  la  believed  on  the  baala  of  experimental  data  that  the  aurfaee 
coefficient  of  heat  transfer,  and  the  adiabatic  wall  temperature  can  be 
approocimated  by  expressions  of  the  form: 

=  0,  t-  k«- 


where  a,  b,  c,  and  d  are  constants 
With  these  two  expressions,  the  boundary  conditions  become: 

•e- «  o  :  “r«  o 

Z  -  O  ;  _  O 

z  «  K  ^»7.  *  ♦A'*-*—  t') 

The  Laplace  transform  operator  Is  now  applied  to  both  the  differ¬ 
ential  equation  and  the  boundary  conditions.  The  operator  is  defined  as 
follows: 

O0 
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The  result  of  this  operatloo  la: 
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M  ^  x^tt.a.  i-bcj 


The  equation  has  thus  been  transformed  to  a  second  order  total  differential 
equation  idiose  most  general  solution  Is: 


T  =  Mp)  ocnJL  %Z 

idiere :  ^ 

The  first  of  the  tao  boundary  conditions  determines  that  B  •  0, 
and  the  second  yields  the  following  first  order  differential  equation: 


— 


a.  (a,st  -h  bcS 
b7’  os«i.  Ck 


6dL 


Or,  writing  everything  In  terms  of  % : 
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c«a.^V. 


rx^c 


^  Ca>^  4*  \» 


^xdL  n 
oT^J 


>  Srirvi 


Noify  let 


AlV)  =  £<.-•'>'■  ft- 1*' *')*»«►<■ 


and  the  differential  equation  becomes: 


H«e  Vi 

_  <1<i]b3 


♦fcc  .  e(a.it  ♦  lid  Z4ft 

b*<i?  t«<*-i*‘ 


or  I  equating  coefficients  of  V“M  C 


dL3^ 

«<■% 


%-  ifytbu.X\K 
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BUdL  4.  bc^  -t.  ^ 


The  homogeaeoue  solution  to  this  equation  can  be  obtained  directly  in  the 
forft  of  an  Integral: 

•►UX4-0K  -I- 

B.h  -  K  ^ 


As  p  approaches  infinity,  the  Laplace  transform  of  the  temperature 
must  approach  zero.  This  can  be  seen  by  the  following  analysis: 


T 


P  m  P  ^ 

Vvr- 


=■  M  +  N 


The  temperature  T  is  bounded  by  the  adiabatic  wall  teiqjerature  at 
all  times.  Hence: 


c«-  + 


and  for 


^  “T^  CL  +  oLL* 


.  40  . 

HonMy  Iqr  tho  theorea  on  upper  bounds  of  integrals: 


N<  i5:i[(»— VysX®'- 


— 

—  If  p  Is  now  allowed  to  approach  inf Inity,  the  products  V.  &  % 

urft*'*^jand  L?  C"^  approach  sero,  and  */p  ,  V^t/^approach  zero. 
Henee  both  M  and  N  approach  zero,  and: 


^  wa  JU^  s  o 

P^OD  F-^0* 


The  hwogeneous  solution  of  however  contains  a  factor  of  ft 

This  factor  produces  a  contribution  to  A  of  the  fora: 


Ah 


Wae 


where  does  not  approach  zero  idien  P  becomes  infinite. 

But  this  contribution  becomes  infinite  when  P  becomes  infinite,  and  henee 
contradicts  the  requirement  that  A  approach  zero.  Thus,  the  coefficient  K 
of  the  homogeneous  solution  must  be  zero.  The  particular  solution  is  found 
by  the  method  of  variation  of  parameters.  A  solution  is  assumed  in  the  fora. 


e’’ 


When  introduced  in  the  differential  equation,  the  function  f  is  expressed 
as  a  quadrature. 
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If  this  expression  is  integrated  by  parts  repeatedly*  an  asymp¬ 
totic  expansion  for  f  is  determined,  which  contains  increasing  negative 
powers  of  ^  •  This  expansion  is  as  follows: 


■¥  {.xyi  4- V\. 


But:  »  I  + 
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Hencei  -JI  •  ^  3lx  k  -  *- 

The  factor  can  non  be  expanded  in  Inverse  powers  of 

,  to  yields  ' 

+  6e.“*^^ -i-\8e“‘^‘^'^  -f  38cr®^''4--^ 
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Introducing  this  into  the  expression  for  f  >  and  arranging 
terms  in  order  of  negative  posers  of  \  ,  the  expression  beecmiest 

“  ^C«^  +  4.01X  4ar)c*'4-]i^ 

+  3 €-*•^4'  4-  3a.  -»•— ^ 

-  h-cibx  4.a.^er^'' 

■f  4-  S’oe-^^’^  4-i7oe"‘'^‘'  4 — 

+  + 2.e“*’^'’  -4- ze"^'*'  +  ■*■•‘^3 

4<-%x4.»)c-*'^  +  ^«x4.i0  e-^""^0z>w+z5-)e‘*i-.] 
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« 

+.  38e“‘*'' -*• ---^1  + --- 

and  making  use  of  the  inverse  Laplace  transforms 

The  following  expression  for  the  temperature  is  derived: 

■T-E^  +  ■'•■"W*-^ I »•’*>•+>■►  “1^1 

'*Mme  J 
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ai»i  where  the  coefficients  C <«•««».  >  are: 

Cea  aa  ^cvsr  S  - 

C|»  a  atiTST  C,^  s-3a.^cot 

Cl  a  a  0.C  VST*  ^ M,A  *  “  "S^ tC 

Cs3  a  acVT*  C,^  a-7a^c•c 

C.a-  -  a.*cn<^  4- iCa.4.-»- Wc^VoT 

C.s-  =  ffo^CoC^  4- -l-ViO'fy  -a.eLbck^ 

Ca.j- =i  V3  a.»Co<>^  4- jlC.«-^ +^<J^'ior-6abcKi^ 

Car  *A5-a.*co^'b  4- ZC<U.  4- bOYSr  -  14-a.Wc.Vi^ 

C04  —  —  ^coc**— Z(a*^  "V  ba^fluot  — 

C/c  »  — Tfl^CBC**  -  +  fac)«bot  -  4-  €9^  6ckK. 

C*0  a  €«<"•- /o<aflt 4- bc)abet  -  *•%  a.4eo{  4*  ao  a^'bc  4.  «C 

C34  a  -M4.‘*C0C*- 30C«.<t4-bC^«L«(— ^4^«-b«0t  4-  30  0,*bc4oC 

C*7  a  a-^C  oc^  4-  Z  (4.a  4-  bO  4-  *Vit  e<^  4-  bdVw' 

C,^  a  3a.^coc'^  4- io(,<bJL4|,c)a.*oc^  4- 4- fibdVST 

-4-(ceL4- be)  bVvV^  —  IJCO.*  btk  -  X  ac  b* K*'V5^ 

C^y  =  -H3.tf(a.d.4>  bC)o.*'P(^a  4-  4  BWdVST 

-  »*.UJ  4kc>bV»-/7'  -  flee?  bcK«<.^*'  -  6cc  b'k^yST' 

C  =  ia.fleifcoc^  4-^oC«.A  +  bc)fl5'ec^*  4.  *^a**bcoc^«'  4-  flbdL-VoT 
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--38Ce*A +bc)bk-flr  -  343  tt?bcK.«4^*  —  1  a.  0.0 b' k^’VK’ 
a  -  O^Ca(^-;4(^4b&)0?K*’‘-  *9k  7(o.^4-bc)^c( 

Co,  »a7co<^  4*(ae4^fcc)<|?e<^  4.  4 

4. /rCo.W.  4- bc)a.bx^  4-ioa*flco^^* 
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■  Q?C  ^Ac)  a^kc  oc^  -w  G  4^  01.^ 

+  3A>Ca.a.  -*.  be)  «*A  oc'^  4.  »•%.  <»?  b*’  CeC*^  +  5*7  0.6*  4.  oC ^ 
4.  aj’CAflt  ♦Ac)  b^oC*^ 

A^«.^  4-  Ac)  o/oc'^  -  o7  Ac  or***-  fi  o/“Ai6  oc» 
-to(*.i  +  fce)o.‘»b«t*-  «».*6*'C.l*-75-*‘4*-a^o<‘ 

—  4*  6c)flL,A**C^— 4-5'O.A*  CoC** 


As  a  check  on  this  result,  xninerlcal  calculations  have  been  made 
for  two  points,  and  compared  with  results  of  a  relaxation  procedure.  The 
series  was  found  to  converge  very  rapidly  with  respect  to  m,  but  not  very 
rapidly  with  respect  to  n.  The  convergence  is  improved  greatly  at  smaller 
values  of  The  calculations  were  made  for  the  following  set  of  condi¬ 
tions; 


Acceleration  •  1  Kach  number  per  30  seconds  starting  from 

Mo  ■  1*4 

OL  r  .909  1/ft 
b  =  136.3  1/ft  -hr 
C  s  26,«00  opAr 
d  =  776,000  ®F/hr2 
Vv.  s  0.20  ft 
0^*  =  .384  ft2Ar 
=  1 


Por  m  .0383  hrs,  or  IS  s  6* 


AT«  T, 


6CO  "F 
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Fop  r  .02165  hra,  or  II  ■  4: 


at  »  \€0.S  •F 


The  relaxation  procedtire  glree  values  of  672  axid  I6d.5^  respectively 
for  these  conditions. 


IJ 

calculated 

AT*relaxatiaga 

%  difference 

H 

160. 5<T 

168. 5°F 

4.75 

[I 

660<^ 

672®F 

1.79 

These  percentage  differences  are  elthin  the  error  of  the  finite 
difference  procedure,  so  that  this  i^ould  serve  as  a  valid  check  on  the 
analytical  procedure. 


p 


48 


APP] 


JLl 


EITECT  OF  DEFLECTIONS  ON  AERODYMAMIC  COEFFICIENTS 

In  this  appendiZy  the  changes  in  aerodynamic  coefficients  due  to 
the  thermal  deflections  are  calculated.  The  work  is  based  upon  the  Ackeret 
thin  wing  approxlaationSf  i.e.,  the  flow  around  every  chordwise  section  of 
the  wing  will  be  considered  as  two-dimensional  and  the  effect  of  the  tip  in 
reducing  the  fluid  pressure  will  be  neglected.  For  the  configuration  shown 
in  figure  11,  the  Ackeret  formulae  for  pressure  coefficients  on  the  upper 
and  lower  surfaces  of  the  wing  are  as  follows: 


c. 


p-p* 


Cm  ^ 


4.  ur  *>•  Wx 

Si*  OteCi-»-\l)  H-UT-  VsL 
But,  or  »  AC\^  4. 


where,  ACX')  •  b.V 

BW  »  4- 


CC1)  -  4-  b,V 


49  - 


Heaoat 


Lift 

'/tVU‘(a4H 


50 


N»ic  0^  Msmcn’t  Mm 


}fc  V  U'  l  a.*')*-  •  I 


T^*  1 5®'*  ■ 


^  5  (f .»  1- 


Dti<i» 


(c»*  -  Cr«.X«,  -*■  '®/i  -*• 


-  51  - 


**-^v4  ■►'•Vi  +  ^*-*/m  ■*■  •♦''•* Vs 


-  52  - 

Pot  s  100  #/ft2 

0(«  s  .0242  radians 

A  8a^}la  oaleulatlon  will  noa  be  nade  for  the  ease  of  a  wing  of 
R  ■  and  OC  s  6.28  x  10*^«  Rdaolta  of  oaleulatlona  for  other  vlnga  are 
tahnlated  in  table  %  ■.• 
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<*'•«  coefficients  sre  defined 
ess  change  In  coefficient  due  to  deflection  ,  ^ 

undeflected  coefficient  * 

%  change  in  Ci,  «  a  Wfff  « |00  «  1.62.% 
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